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ABSTRACT 

We present a compact parametric representation of the smooth bright muhisolution 
solutions for the modified Camassa-Holm (mCH) equation with cubic nonlinearity. We 
first transform the mCH equation to an associated mCH equation through a reciprocal 
transformation and then find a novel Backlund transformation between solutions of the 
associated mCH equation and a model equation for shallow-water waves (SWW) intro- 
duced by Ablowitz at al. We combine this result with the expressions of the multisoliton 
solutions for the SWW and modified Korteweg-de Vries equations to obtain the multi- 
soliton solutions of the mCH equation. Subsequently, we investigate the properties of the 
one- and two-soliton solutions as well as the general multisoliton solutions. We show that 
the smoothness of the solutions is assured only if the amplitude parameters of solitons 
satisfy certain conditions. We also find that at a critical value of the parameter beyond 
which the solution becomes singular, the soliton solution exhibits a different feature from 
that of the peakon solution of the CH equation. Then, by performing an asymptotic 
analysis for large time, we obtain the formula for the phase shift and confirm the solitonic 
nature of the multisoliton solutions. Last, we use the Backlund transformation to derive 
an infinite number of conservation laws of the mCH equation. 
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I. INTRODUCTION 

A modified Camassa-Holm (mCH) equation with cubic nonlinearity is a current inter- 
est in the theory of sohtons. It may be written in the form 

mt + [m{u^ - ul)]x ^ 0, m^u-u^^, (1.1) 

where u — u{x, t) is a real- valued function of time t and a spatial variable x, and the sub- 
scripts X and t appended to m and u denote partial differentiation. The mCH equation 
was found independently by several researchers using a procedure that generates new in- 
tegrable systems from known integrable bi-hamiltonian systems. It is well-known that 
the method yields the CH equation^ when apphed to the Korteweg-de Vries (KdV) equa- 
tion whereas it yields Eq. (1.1) when applied to the modified KdV (mKdV) equation. -"^"^ 
It also arises from the physical system like nonlinear water waves. ^'^ Although a large 
number of works have been devoted to the CH equation,® the study of the mCH equation 
is still under way. Specifically, a recent work by Qiao reveals that it admits W/M-shape 
cusp solitons unlike the peakon solutions of the CH equation. ^'^ Quite recently, smooth 
dark soliton solutions of the integrable hierarchy including the mCH equation were ob- 
tained for the variable m using the inverse scattering transform (1ST) method and the 
properties of the one- and two-soliton solutions were explored.^ 

The purpose of the present paper is to construct the bright A'"-soliton solutions of 
the mCH equation under the boundary condition — > as |x| — )■ oo, where N is an 
arbitrary positive integer and uq is a positive constant. We employ an exact method of 
solution which worked effectively for constructing the A^-soliton solutions of the CH and 
Degaperis-Procesi equations. ^^^^ As a byproduct, an infinite number of conservation laws 
are produced in a simple manner without recourse to the 1ST. 

The paper is organized as follows. In Sec. II, the mCH is transformed to an asso- 
ciated mCH equation by introducing a reciprocal transformation. We find a remarkable 
Backlund transformation connecting solutions of the associated CH equation with those 
of a model equation for shallow-water waves (SWW equation for short) introduced by 
Ablowitz et al}^ This allows us to obtain the parametric representation for the bright 
A'"-sohton solution of the mCH equation when coupled with the A^-sohton solutions of 
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the SWW and mKdV equations. In Sec. Ill, the properties of the one- and two sohton 
solutions are investigated in detail. In particular, we obtain a condition necessary for the 
existence of smooth solitons and show that there exists a critical value of the amplitude 
parameter of the soliton beyond which the solution becomes singular. Subsequently, the 
asymptotic behavior of the A^-soliton solution is briefly discussed, whereby we derive the 
formula for the phase shift and confirm the solitonic behavior of the solution. In Sec. IV, 
an infinite number of conservation laws are derived simply through the Backlund trans- 
formation and their explicit expressions are presented for the first few members. Section 
V is devoted to concluding remarks. 

II. CONSTRUCTION OF MULTISOLITON SOLUTIONS 
A. Associated mCH equation 

We first introduce a coordinate transformation {x, t) — >■ {y, r) by 

dy — mdx — m{u^ — ul)dt, dr — dt, (2.1a) 

subjected to the restriction m > 0. Accordingly, the x and t derivatives can be rewritten 
as 

d d d d / 2 2\ "9 /o 7\ 

d^^'^w dt=d^-'^^'' -""^W ^ ^ ^^ 

It then follows from (2.1) that the variable x — x{y, r) satisfies a system of linear partial 
differential equations (PDEs) 

Xy — — , Xr — — rn^ul. (2.2) 
m " 

We apply the transformation (2.1) to the mCH equation and find that it can be recast 
into the form 

+ 2m^Uy = 0. (2.3) 

The integr ability condition x-^y — Xyr of the system (2.2) is satisfied automatically by 
virtue of Eq. (2.3). Using (2.3) and the relation u^x — rri^Uyy -\- mmyUy, the variable 
u{— m + Uxx) is expressed in terms of m as 

u^m+l;m(^] . (2.4) 

ry 
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Last, if we define the new variable p by p — 1/m and substitute u from (2.4) into (2.3), 
we obtain a nonlinear PDE for p: 



PPryy " PyPry " pV " '^Py = 0. 



(2.5) 



This evolution equation for p in the independent variables r and y is the reciprocal 
transformation for the mCH equation. We call it the associated mCH equation. 

B. Multisoliton solutions 

The main result in this paper is given by the following theorem: 
Theorem 2.1: The mCH equation (1-1) admits the parametric representation for the 
N-soliton solution 

U = U{y, t) =Uo- (In /7)rjy, 



X = x(y, r) = — + uIt + 2 In^ + d, 
uo J 



where f and f are tau-functions given by 



f = ^ exp 

1^=0,1 



N 



l<i<j<N 



i=l 



i"=0,l 



N 



l<i<j<N 



i=l 



with 



y- 



(uoki) 



:T — 



Un 



- Akf \ s-yio 



l,2,...,iV), 



1 - Uoki 



1 + Uokj 



, (i = l,2,...,Ar), 



QYij — 



Jv-i k'i 



, {i,j^l,2,...,N;i^j). 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



(2.10a) 

(2.106) 
(2.10c) 



~\~ kj 

Here, ki andyio are the amplitude and phase parameters of the ith soliton, respectively, 
d is an arbitrary constant, s is an auxiliary time variable and the notation X^^=o i if^P^'i^s 
the summation over all possible combination of /ii — 0, 1, 112 — 0, 1, /in — 0, 1. 
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The parametric solution provided by Theorem 2.1 would give rise to the smooth bright 
A^-soliton solution on a constant background u — uoif the conditions < uoki < V3/2 {i — 
1, 2, N) are imposed on the amphtude parameters. The variable s introduced in (2.10a) 
plays the role of the time variable for both the KdV and mKdV equations. It can be set 

to zero after all the computations have been completed. 

Now, Theorem 2.1 follows straightforwardly from Propositions 2.1 and 2.2 below when 
combined with the system (2.2), which we shall now demonstrate. 

Proposition 2.1: The variables p., q and r defined by 

p=- + 2(\J-\ , 

«0 V fjy 

g = -2 (In f)yy, 

= 4 - 8(ln/)2,j;, 

satisfy the nonlinear PDEs 

Ps + Gp'^Py - Apyyy = 0, 

Qr + 2ul Qy + qqr + 2ul qydy^qr - ul qryy = 0, {d'^ ^ - dy), 

+ Q'Ty - A'lyyy = 0, 

respectively. Then, the variables r and q are expressed in terms of p as 

r^p^ + 2py, 

^ 2^4 4x^2- 

Proof: Equation (2.14) is the mKdV equation of defocusing type. It exhibits the dark 
A'"-sohton solution (2.11) with the tau-functions (2.8) and (2.9). -^^'-^^ Equation (2.15) is a 
SWW equation for which the A'"-soliton solution is given by (2.12) with the tau-function 
(2.8) whereas Eq. (2.16) is the KdV equation whose tau-function has the same functional 
form as that of the SWW equation except the time dependence. ^^'^^ The latter follows 
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(2.11) 

(2.12) 
(2.13) 

(2.14) 
(2.15) 
(2.16) 

(2.17) 
(2.18) 



from the fact that the SWW equation belongs to a member of the integrable hierarchy of 
the KdV equation. 

The relation (2.17) is the Miura transformation^^ which connects solutions of the 
mKdV equation with those of the KdV equation. Actually, a direct substitution of (2.17) 
into the left-hand side of (2.16) reveals that 

r, + 6rry - Aryyy = 2 + {p, + Gp'^Py - 4:pyyy). (2.19) 

The relation (2.18) stems simply from (2.12), (2.13) and (2.17) by eliminating the 
tau-function /. □ 

The relation (2.18) is a Backlund transformation between solutions of the associated 
mCH equation and the SWW equation, as will be shown by Proposition 2.2. If one 
substitutes (2.11)-(2.13) into (2.17) and (2.18), then one can see that they reduce to the 
following bilinear equation for / and /': 

4/ - 2/;/. + ffyy + ^ifyf - ffv) = 0. (2.20) 

The above equation is a constituent of the system of bilinear equations for the mKdV 
equation and hence the tau-functions / and /' from (2.8) and (2.9) solve Eq. (2.20).^^ 

Proposition 2.2: Let 

Q^q^ + 2ul Qy + Aul qqr + 2ul Qyd'^qr - ul Qryy, (2.21) 

j^^Pryy_ PyP^ - pp, - 2 ^. (2.22) 
p p^ p^ 

Then, Q and R are connected by the relation 

under the transformation (2.18). 

Proof: Substituting (2.18) into (2.21), one can recast Q into the form 

Q = Y^iPyy + PPy){'^Uo + dy^ppr)} + p'^Pry + ^PPrPy + P^'Pr - Pryyy - PPryy]- 



On the other hand, integration of (2.22) with respect to y under the boundary condition 
p — >■ 1/uo, \y\ — >■ oo gives 

2uo + d-\ppr) = -d-^R + ^{pry + 2), 

whereas a direct computation using (2.22) leads to 

(v^R) 1 

— ^ +p'^R= -p^Pry - SpPrPy - P^Pr + Pryyy + PPryy - - {Pyy + PPy) {Pry + 2) . 

The relation (2.22) follows by introducing above two expressions into Q. Actually, 

Q = -^[iPyy+PPy)d;'R+pRy + {2py+p')R] = [d^ + P) d^^P^v'^^' 
This completes the proof. □ 

It is obvious from (2.23) that the equation R — yields the equation Q — 0. But, the 
converse statement is not true in general. For the A^-soliton solution (2.11) which satisfies 
the boundary condition p po (= I/uq), \y\ — )■ oo, however, p is a solution of Eq. (2.5), 
or i? = if g from (2.12) satisfies Eq. (2.15), or Q = 0. To see this, assume first that 
Q = 0. Integrating the resultant equation from (2.23) with respect to y and using the 
boundary condition p ^ pQ,y ^ +oo, one obtains 

d 



-^{pdy^R) = Ro{r)exp 



/■oo 

/ {p-Po)dy-poy 



(2.24) 



where Ro{t) is an integration constant. A slight modification of the relation for dy ^R 
derived in the proof of Proposition 2.2 gives 

o 1 ^ Ptv 1 5 Z"^/ 2 2x , 2 2 

" p 2dT Jy P Po 

As confirmed easily by integrating (2.18) with respect to y and using (2.12), the integral 
X^(p^~Po)^y becomes a constant independent of r. It turns out that limj,_^_oo dy^R — 0. 
With this fact in mind, the left-hand side of (2.24) vanishes in the hmit of y — > — oo. On 
the other hand, the right-hand side diverges in this limit. Thus, one must put Rq — 0, 
which yields {pdy^R)y = 0. Integrating again with respect to y, one has dy^R = Ri{t)/p, 
where Ri{t) is an integration constant. In view of the boundary condition, this constant 



must be zero, and hence dy^R — 0. Differentiating this relation by y, one finally arrives 
at the relation R — 0. 



We are now ready for proving Theorem 1.1. 

Proof: First, referring to (2.2) and (2.11), one obtains 



Integrating this relation with respect to y leads to 

uo J 

where (i(r) is an integration constant. Differentiate the above expression by r and use 
(2.2) to find the relation 



This expression reduces, upon taking the limit y — )■ oo and substituting the limiting 
values u — >■ uo,Uy — >■ 0, [ln(/Y/)]j^ 0, to d'{T) — Kq, which gives d{r) — u^t + d after 
integrating with respect to r. This proves (2.7). 

It now follows by integrating Eq. (2.3) with respect to y and substituting from 
(2.18) into the resultant expression that 



III. PROPERTIES OF SOLUTIONS 

A. One-soliton solution 

The tau-functions / and /' corresponding to the one-soliton solution are given by 
(2.8), (2.9) and (2.10) with = 1. They read 





Last, substitution of (2.11) and (2.12) into the above expression yields 




which is just (2.6). 



□ 



/ = l + e^ 



(3.1) 
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/' = l + e«-^ (3.2) 

with 

P = hill — CT — iirX r = 

1 - {uoky 



^ = k{y-cT-yo), c=- — , (3.3a) 



e-^ = (3.36) 

where we have put ^ — ^i,(p — (f)i,k — ki and yo — yio for simphcity and set s — 0. We 
assume A; > hereafter. 

The parametric representation of the one-sohton solution foUows by introducing (3.1) 
and (3.2) into (2.6) and (2.7). It may be written in the form 

k'^c{acosh{^-^o) + 1} 
{cosh(^ - ^o) + ar 

X — ct — xq = — - + 2 In , — , , (3.46) 

where 

1^0 1 - [uQky 

and we have put d — XQ — yQ/ UQ—\n[{l — a) / {l-\-a)\. Note that cfrom (3.4d) isthe vefocity 
of the sohton in the (x, t) coordinate system. The smoothness of the sohton solution is 
assured by the condition m = 1/p > (see (2.1)) which imposes certain restriction on the 
parameter a. To show this explicitly, we compute Xy[= p) from (3.4b) and obtain 

1 2ka 
P— , . (3.5) 

Uo y/T^^ COsh(e - 6) + 1 

Recall that p is a dark soliton solution of the mKdV equation (2.14). The required 



condition is then found to be satisfied if I/uq — 2ka/{-\/l — + 1) > 0. Thus, the 
parameter a must lie in the interval 

< a < (3.6) 

One can see from (3.4) that the one-sohton solution represents a bright soliton on 
a constant background u = Uq whose center position Xc is located at Xc = ct + Xq + 
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-6 -4 -2 2 4 6 

X 

FIG. 1. The profile of smooth sohtons. a = 0.5 (dashed curve), a = 0.7 (dotted curve), 
a = \/3/2 (sohd curve). 

^o/iuok) (^ = ^o)- In view of this observation, the amphtude of the sohton with respect 

to the background field, which we denote by A, amounts to 

A^-t^il-Vl^'). (3.7) 

V 1 — « 

Eliminating the parameter a from (3.4d) and (3.7) leads to the amplitude- velocity relation 
of the soliton 

A = ^2{c-ul) - 2uo. (3.8) 
The inequality (3.6) restricts possible values of c and A. Explicitly, 

3ul<c< 9ul, 0<A< 2uo. (3.9) 

If the condition (3.6) is satisfied, then the one-soliton solution represents a smooth 
bright solitary wave travelling to the right with velocity c and amplitude A. As a increases, 
the amplitude grows and the width narrows. This observation is illustrated in Fig. 1, 
whereby the profile of smooth solitons with a background = 1 is depicted agianst the 
stationary variable X — x — ct — xq ior three representative values of a. Of particular 
interest is the limiting profile of the soliton when the parameter a tends to the upper 
hmit a = \/3/2 (~ 0.866) of the inequality (3.6) (see sohd curve in Fig. 1). In this limit, 
the smoothness of the profile would be lost at the center position (or the crest) of the 
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-6 -4 -2 



FIG. 2. The profile of singular solitons. a = 0.94 (dashed curve), a = 0.96 (dotted curve), 
a = 0.97 (solid curve). 



soliton. The typical profile is also depicted in Fig. 2 for singular solitons which violate 
the condition (3.6). 

To explore the feature of the solution in more detail, we shall investigate the profile 
of the soliton near its crest. To this end, we first derive the approximate expression of 
X from (3.4b). If we choose the phase parameter of the soliton appropriately such that 
^ = ^0 corresponds to X — 0, then the power series expansion of X with respect to ^ — 
is found to be as 



X = — (1-2^1^) (e -Co) + ^^^^(l-v^r^) (C - Co)' + o ((e - eo)') . 



In particular, for a — -\/3/2, (3.10) reduces to 



^ = ^ (e - 6)' + O ((e - Co)') 



(3.10) 



(3.11) 



On the other hand, the expansion of u from (3.4a) reads 



u — Uq 



1 + 



1 + a 
which simplifies to 



6 3(-2 + a),, ^ 16-13a + a2 
+ ). . J (e - 6)' + 



{1 + ay 



4(1 + a)^ 



(e - 6)' + O ((e - ^of) 



U — Uq 



3--(e-6r + 0((e-eo)') 



(3.12) 
(3.13) 
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when a — 2{a — -\/3/2). Note that the coefficient of — ^o)^ in the expansion (3.12) 
vanishes for this specific value of a. By inserting ^ — from (3.11) into (3.13), we find 
that the profile of the soliton near the crest is approximated by the expression 



3 _ ^ X'/' 
2 



(3.14) 



Thus, we can see that the smoothness of the solution is lost at the crest X = 0. More 
precisely, u has only a continuous first derivative near the crest and the nth derivative 
does not exists for n > 2 at the crest. This novel feature is in striking contrast to the 
appearance of the singularity (or discontinuous first derivative at the crest) in the peakon 
solution which can be deduced from the smooth soliton solution of the CH equation 
through an appropriate limiting process. ^^'^^ 

B. Two-soliton solution 

The tau-functions / and /' for the two-soliton solution can be written as 



/ = 1 + e^i + e^2 + e^^+^^ (3.15a) 

Vfci + k2j 

\ki + k2j 



where 



2u 



3 



= ki{y -CiT- Vio), Ci = ^ _ (^^°^.)2 ' = 1> 2), (3.15c) 

e-^^ = ^^^, (^ = l,2). (3.15ci) 
1 + uoki 

We investigate the asymptotic behavior of the two-sohton solution for large time and 
show that it describes the elastic interaction of two bright sohtons on a background field 

u = uq. To this end, let q (i = 1,2) be the velocity of the ithe soliton in the {x,t) 
coordinate system and assume that < C2 < Ci. In addition, we impose the conditions 
< Uoki < 73/2, {i = 1,2) to assure the smoothness of the solution. 

First, we take the limit t — )■ — cxo with being fixed. In this limit, ^2 —00. Then, 
the leading-order asymptotics of the tau-functions are given by 

/~l-Fe^S (3.16a) 
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It folows from (2.6) and (2.7) and (3.16) that 

kfciiai cosh({i - ^lo) + 1} 



{cosh(^i - ^lo) + aiY 



X — Cit — Xio ~ — ; — him 



where 



1 



- i In I ilL^^ 



^ ^ C]_ _^ ^2 ^ -»o{3 - {upkif} 
Uq ° 1 - {uoki)'^ 



(3.166) 

(3.17a) 
(3.176) 

(3.17c) 
(3.17c^) 



In the limit t — )■ +oo, on the other hand, ^2 — > +00. The expressions corresponding 
to (3.16) and (3.17) read 



€2— <^2 



1 + 



ki - k2 



X — Cit — XiQ 



ki + k2^ 

klcijai cosh(ei - 60 + Sj^^) + 1} 
{cosh(a-eio + 5^^) + ai}2 ' 

Uoh ^ I y/lT^i e^i -«io+4+^ + y/T^^ J 



where 



In 



ki - k2 

ki + k2 



Let be the center position of the ith soliton. Then, as i — >^ —00, we find 

Xic Cit + Xio + (6=Clo)- 

As t — >■ +00, on the other hand, x\c reads 

Xi, ~ Cii + Xio + ^(60 - ^t^) - 202, (6 = Cio - 



(3.18a) 

(3.186) 
(3.19a) 
(3.196) 

(3.19c) 
(3.20) 



(3.21) 



13 



3 
2.5 

2 

=1 1.5 
1 

0.5 




100 



-50 



50 



100 



FIG. 3. The interaction of two sohtons. t = —10 (left: before collision), t = (middle: 
during collision), t = 10 (right: after collision). 

The above analysis shows that the asymptotic state of the solution for large time is 
represented by a superposition of two single solitons in the rest frame of reference. The net 
effect of the interaction between solitons is the phase shift, which we shall now evaluate. 
To this end, we define the phase shift of the ith soliton by 

Ai ^ Xic{t ^ +oo) - Xic{t ^ -oo), (i = l,2). (3.22) 

Then, we see from (3.20) and (3.21) that the large soliton suffers a phase shift 

Ai = -^In r ^^'| ' - In f l±^'l ' . (3.23) 



Uoh \k1-\-k2J \I-U0k2j 
Performing the similar asymptotic analysis for the small soliton, we obtain the formula 
for the phase shift. We quote only the final result. 



1 , fki-h 
-In 



1 + 



+ ln . (3.24) 

Uok2 \ki + k2j V 1 - uoki J 

The first term of the above formulas is the same as the phase shift arising from the 

interaction of two solitons for the KdV equation and the second term comes from the 

mapping (2.7). A new feature appears due to this additional term. An inspection reveals 

that under the condition C2 < ci {k2 < ki), Ai is always positive whereas A2 takes any 

sign depending on values of ki and /c2. Specifically, there exists a critical curve along which 
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Ai = A2 and beyond which Ai < A2 so that the small soliton exhibits a larger phase 
shift than that of the large sohton. It is interesting to observe that the formulas (3.23) 
and (3.24) coincide formally with those of the two-sohton solution of the CH equation if 
one identifies Uq with the parameter k which enters into a linear dispersive term 2k^Ux in 
the CH equation.^ 

The interaction process of two solitons is illustrated in Fig. 3 for three values of time. 
In this example, Uq = l,Ci = 0.8, C2 = 0.5, = 1.74, A2 = 0.83. The phase shifts of 
the large and small solitons evaluated from (3.23) and (3.24) amount to Ai = 2.20 and 
A2 = —1.89, respectively. 

C. A'-soIiton solution 

We discuss briefiy the asymptotic behavior of the general A'"-soliton solution. First, 
we order the magnitude of the velocity of each soliton as < < c^-i < ... < Ci and 
impose the conditions < Uoki < \/3/2, (i = 1, 2, A'"). We take the limit — >■ —00 
with the coordinate of the ith soliton being fixed. Then, ^1, {2, 6-1 ~^ +00, and 
Ci+i,C,i+2, ■■■jCn —00. The asymptotic forms of / and /' from (2.8) and (2.9) are found 
to be as 




(3.25a) 




2 



exp 5^(0-0,) (e^-'^^+^'"') , 



(3.256) 



Lj=i 



where 




(3.25c) 



Substitution of (3.25) into (2.6) and (2.7) yields 



Ci{ai cosh(^j - + Sl + 1} 



(3.26a) 



X Cjt XjQ 



+ 2 In 



a/1 - a^e' 



2Y,4>r (3.266) 



i-1 



sjl + ai e' 
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In the limit — > +00, on the other hand, ^1, ^2, Ci-i ~^ —00, and ^j+i, 
+00 and the expressions corresponding to (2.25) and (2.26) are given by 



/' 



n 



n 



kj ki 



i+l<j<l<N ^' ^ 
2 



kj + ki 



exp 



N 



U=i+1 



(^l + e^M^'), (3.27a) 



i+l<j<l<N 



kj ki 
kj + ki 



exp 



N 



( 



1 + e^' 



{+) 



where 



and 



j=i+i 



ki ~\~ kj 



kfc,{ai cosh(^i - + 5, 



(+)^ 



1} 



{cosh(^i - 



iO 



5. 



+ aiy 



^^ ,01/ ^r^e«>-««+'^>*^' + ^TT a.; 

X — Cit — Xio ~ — ; — h 2 in 



We can see that the asymptotic form of the A^-sohton solution is a superposition of 
N single sohtons each of which has the form given by (3.4). The phase shift of the ith 
soliton can be derived from (3.26b) and (3.28b). It reads 

■j-i 



TV 

2 E 

j=i+i 



{3.27b) 



(3.27c) 



(3.28a) 



(3.286) 



1 

Uoki 



ki kj ^ f ki kj 



ki ~\~ kj 



N 



j=i+l 
2 



ki ~\~ kj 



j = l ^ ^ ■> ^ j=i+l ^ 



Uok^ 



(3.29) 



In the special case of iV = 2, the formulas (3.29) reduce to (3.23) and (3.24). They clearly 
show that each soliton has pairwise interactions with other solitons, i.e., there are no 
many-particle collisions among solitons. This feature is common to that of the A'^-soliton 
solutions of integrable nonhnear PDEs. 

IV. CONSERVATION LAWS 



The Backlund transformation (2.18) between solutions p and q of the two integrable 
equations allows us to construct an infinite number of conservation laws of the mCH 
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equation in a simple manner. To demonstrate this, we first note that the SWW equation 
(2.15) exhibits local conservation laws of the form 

Wn,T^jn,y, (^ = 0,1,2,..), (4.1) 

where the conserved density Wn and associated flux j„ are polynomials of q and its y- 
derivativcs. Wc rewrite (4.1) in terms of the variables x and t by using (2.1). Substituting 
Eq. (1.1) into the resultant equation, we obtain 

{mwn)t = [jn - m{u^ - ul)wn]x- (4.2) 
It turns out that the quantities 

/oo 
mWndx, (n = 0, 1,2, ...), (4.3) 
-oo 

become the conservation laws of the mCH equation if one introduces q from (2.18) into 

Wn- 

The conservation laws of the SWW equation take the same form as those of the KdV 
equation since both equations belong to a common integrable hierarchy. For complete- 
ness, we reproduce them shortly. We start from the Lax pair associated with the SWW 
equation^° 

i^yy - # = — ^, (4.4a) 
(4A - l)^Pr = 2mo(«o + d-^qr)^y - ulqriJ, (4.46) 
where A is a spectral parameter. If we put w — ipy/ip, we can rewrite (4.4) into the form 

Wy + w^ — q=^, (4.5a) 

Uq 

(4A - l)wr = [2ul{uo + d-^qr)w - ulqr]y. (4.56) 

Thus, the quantity w dy is conserved in r. Expanding w in powers of e(= Uo/{2VX)) 
ets w — l/(2e) + Yl'^=i ^"""^n and substituting it into (4.5a), we obtain, after equating the 
coefficients of e", the recursion relation that determines Wn 

n-l 

Wn+1 = Wn,y - ^ Wn-mWm, {n>2), (4.6) 
m=l 
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with the initial conditions Wi — q and W2 — —Qy We quote a few conserved densities: 



^^3 = Qyy - g^ W4 = {-Qyy + 2q\, W5 = {Qyyy " 6qqy)y + 2q^ + qy. (4.7) 



The quantity q from (2.18) can be rewritten in terms of m as 



+ 



(4.8) 



2m^ Aw? 

Last, substituting (4.7) with (4.8) into (4.3) and taking into account the relation 
d/dy = m'^d/dx, we obtain a sequence of conservation laws of the mCH equation ex- 
pressed in terms of the original spatial variable. We see that Wn for even n takes the form 
of a perfect derivative {Wn)y and hence the quantity /„ for even n vanishes identically. 
Actually, J^^m{Wn)ydx — J^^{Wn)xdx — 0. Hence, only the quantities /„ for odd n 
survive. We can express the nontrivial conservation laws as hn+i — X]m=o ^nmlrm where 
Unrn ^rc coustauts depending on uq. The first four of read as follows: 



J — c 



(m — uo)dx, 



) dx, 



1 , mi. 
1_ 4 — E. 



Un 



dx, 



1 + 8^ + 20^ 



56 4 



dx. 



(4.9a) 

(4.96) 
(4.9c) 
{AM) 



Uq m' m' 
Note that Jq follows directly from the mCH equation (1.1). 

Last, we remark that an infinite number of conservation laws have also been obtained 
quite recently using the fact that the mCH equation describes pseudo-spherical surfaces^^ 
which is equivalent to the existence of the Lax representation for the mCH equation.^ On 
the other hand, our derivation is based on a purely algebraic procedure relying on the 
Backlund transformation. 



V. CONCLUDING REMARKS 



In this paper, we have developed a systematic procedure for constructing smooth 
soliton solutions of the mCH equation. In the process, the associated mCH equation has 
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played a central role which has been derived from the mCH equation via a reciprocal 
transformation. It is important that certain condition must be imposed on the soliton 
parameters to assure the smoothness of the solutions. Specifically, a detailed inspection of 
the one-soliton solution reveals that there exists a critical value of the amplitude parameter 
beyond which the solution becomes a many- valued function. This feature differs from that 
of the onc-soliton solution of the CH equation for which the pcakon solution is produced 
from the smooth soliton solution in the limit of small dispersion. ^^'^^ Whether the singular 
solutions such as the W/M-shape soliton and cusp soliton found by Qiao^ can be obtained 
from smooth solitons by means of appropriate limiting procedure is an interesting issue 
to be studied in a separate context. 

An exact method of solution presented here is also applicable to a variant of the mCH 
equation, for example, 

mt + [m{u^ - ul)]x + = 0, m^u- u^x, (5.1) 

where 7 is a real parameter. When 7 = 0, this equation is shown to exhibit peakons 
under the boundary condition u ^ as \x\ ^ 00.^^ In addition, another version of the 
mCH equation with cubic nonlinearity 

rrit + u^nix + 2>uUxm = 0, m = u — Uxx, (5.2) 

attracts an interest because of its different mathematical structure from that of the mCH 
equation under consideration. This equation has been discovered by Novikov^^ and admits 
a Lax pair related to a negative flow of the Sawada-Kotera hierarchy through a reciprocal 
transformation.^"^ For both equations (5.1) and (5 .2), we obtained smooth single soliton 
solutions which satisfy the boundary condition u ^ uq ets \x\ ^ 00. However, the general 
A^-soliton problem still remains open and it reserves a future study. 

In conclusion, we comment on a paper by Ivanov and Lyons. ^ Applying the 1ST to an 
initial value problem of the mCH equation under the boundary condition m — > mo(= uq) 
as I a; I — )■ 00, they presented the explicit parametric representations of the one- and two- 
soliton solutions for the quantity m(= u — Uxx) which take the form of dark solitons. 
See Fig. 1 in their paper, for example. However, the expression (37) for x is incorrect. 
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Actually, the sign of the first term on the right-hand side of (37) should be minus in place 
of plus. This discrepancy has been caused by an error in Eq. (32). Specifically, a factor 
Qiky should read e""**^^. The same error happened in Eq. (37) as well for the two-soliton 
solution. If these errors were corrected appropriately, then their result would coincide 
with that given in the present paper. Note, however that Ivanov and Lyons obtained the 
expression of m only. On the other hand, we have a compact parametric representation 
of u for the first time. 
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